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The  information  on  this  page  is  offered  to  assist  in  interpreting 
Figures  1-4. 


Minimum 

Principal  Stress/Strain 


Numerical  Value  of  Average 
Stress/Strain 

__  Maximum 

Principal  Stress/Strain 


Mohr's  Circle 

Circle  diameter  is  proportional  to 
the  maximum  shear  stress/ strain 


Twice  the  angle  between 
the  maximum  principal 
stress/strain  and  the 
x-axis 


Notes: 


1.  At  some  x,y  no  stress  circle  is  found  because  the  computational 

procedure  for  stress  did  not  converge.  This  occurred  principally 
along  boundaries  where  the  strain  calculation  may  be  less  accurate. 


2. 


Principal  stresses  are  assumed  to  have  the  same  direction  as 
principal  strains. 


FATIGUE  CRACK  TIP  STRESS  RANGE  AND  ENERGY  CALCULATIONS 
AS  DERIVED  FROM  STRAINS  BY  STEREOIMAGING 


I.  INTRODUCTION 

The  computation  of  stresses  from  the  strains  measured  near  fatigue 
crack  tips  by  stereoimaging  has  been  accomplished.  This  report  outlines 
the  methods  used  to  make  these  calculations,  together  with  the  assumptions 
implicit  in  them,  and  presents  the  results  which  have  been  obtained  to 
date. 


Another  method  has  been  devised  to  calculate  the  energy  expended 
in  creation  of  a  unit  area  of  new  crack,  d^/dA.  The  plastic  strain  range, 
as  determined  by  stereoimaging,  is  used  together  with  the  hysteretic 
energy  loss  per  cycle  to  make  this  calculation.  The  hysteretic  energy 
loss  per  cycle  has  been  both  directly  measured,  and  also  derived,  using 
previously  obtained  subgrain  size  information,  giving  equivalent  results. 

The  information  on  the  effect  of  water  vapor  on  fatigue  crack 
plasticity  which  has  been  determined  during  this  past  year  and  previous 
years  has  been  used  to  include  both  environmental  effects  and  threshold 
effects  into  the  damage  accumulation  model  for  fatigue  crack  growth. 
Finally,  this  report  will  summarize  the  effects  of  water  vapor  on  fatigue 
crack  tip  stresses,  strains,  energies,  and  opening  modes,  and  will  offer 
some  possible  mechanistic  interpretations  of  these  results. 


II.  DETERMINATION  OF  CYCLIC  STRESS  RANGE 
A.  Introduction 


The  cyclic  stress  range  distribution  within  the  fatigue  crack 
tip  plastic  zone  was  previously  determined (D  using  subgrains  as  small 
"stress  gauges."  Because  of  the  successful  correlation  between  strain 
distribution,  as  determined  from  subgrains,  with  that  determined  by  stereo¬ 
imaging  analysis, (2)  it  is  desirable  to  determine  stress  from  stereoimaging 
strain  in  order  to  make  a  similar  correlation.  Another  reason  for  computa¬ 
tion  of  stresses  is  to  try  to  determine  experimentally  the  magnitude  of  the 
mean  stress,  or  hydrostatic  stress,  within  the  plastic  zone.  This  is 
useful  because  of  the  role  mean  stress  is  thought  to  have  in  environmentally 
assisted  crack  propagation. 

B.  Calculation  of  Principal  Stress 

The  surface  of  any  specimen  is,  by  definition,  in  a  state  of 
plane  stress,  which  means  no  constraint  exists  in  the  through-thickness 
direction  which  can  support  a  stress  perpendicular  to  the  specimen  surface 
plane.  It  is,  therefore,  not  possible  to  determine  stresses  in  the  con¬ 
dition  of  plane  strain,  although  all  evidence  collected  to  date  Indicates 


that  the  plastic  zone  sizes  and  shapes  are  about  what  would  be  expected 
for  the  crack  tip  plastic  zone  to  be  in  plane  strain.  Thus,  a  peculiar 
set  of  conditions  pertains:  the  material  constraints  are  such  that  the 
crack  tip  plasticity  is  indicative  of  conditions  within  the  specimen,  so 
that  the  through-thickness  variations  are  minor,  while  the  stresses  must 
be  considered  as  though  the  specimen  itself  is  wholly  in  a  state  of  plane 
stress.  The  significance  of  these  factors  and  their  exact  meaning  are 
still  being  pondered.  However,  our  tentative  conclusion  is  that  for  most 
of  the  plastic  zone,  the  stresses  as  well  as  strains,  reflect  conditions 
on  the  interior  of  the  specimen  as  well  as  on  the  surface;  which  is  to 
say  that  for  the  constrained  plasticity  conditions  used  here,  the  computed 
stresses  are,  more-or-less,  those  which  would  be  found  in  the  specimen 
interior,  although  immediately  at  the  crack  tip,  this  situation  has  a 
larger  degree  of  uncertainty. 

The  strains  measured  by  stereo imaging  are  considered  to  be 
total  strains;  thus  ATi,  the  measured  principal  strain  range,  is 

AT^  “  Ae^  +  Ae^  (1) 

where  Aei  ■  the  plastic  strain  range  and  Ae^  *  the  elastic  strain  range. 
With  several  assumptions.  Appendix  A,  it  is  possible  to  separate  elastic 
and  plastic  strain  range,  and,  through  the  use  of  the  cyclic  stress-strain 
curve,  to  compute  the  magnitude  of  the  principal  stress  ranges  from  the 
principal  strain  ranges.  It  is  also  assumed  that  the  principal  stresses 
have  the  same  direction  as  the  principal  strains. 

Computations  of  the  stresses  for  specimens  cycled  in  the  wet 
and  dry  environments  are  shown  in  Figures  1  and  2.  These  are  Mohr's  circle 
diagrams  of  stress,  which  match  the  Mohr's  circle  diagrams  of  strain,  shown 
in  Figures  3  and  4  (reprinted  from  last  year's  report).  The  angular  varia¬ 
tion  in  effective  stress  range,  Aaeff,  is  shown  for  each  environment  in 
Figures  5  and  6  (see  Eq.  A7  in  Appendix  A  for  definition  of  Aoef f ) .  Also 
shown  on  these  figures  is  the  angular  variation  in  stress  range  computed 
from  subgrain  formation.  There  is  remarkable  agreement.  Table  1,  in  the 
maximum  value  of  stress  at  the  crack  tip  as  determined  by  the  two  different 
methods,  for  both  environments,  even  though  there  are  differences  in  the 
angular  distributions  of  stress.  These  differences  have  not  yet  been 
resolved,  but  they  may  stem  from  changes  in  the  plane  stress-plane  strain 
mix  measured  by  the  two  techniques.  (Stereoimaging  must  be  done  on  the 
specimen  surface,  while  the  subgrain  data  is  best  taken  from  the  near 
subsurface  region;  a  depth  fr,om  the  surface  of  0.1  mm  was  used  in  previous 
work.)  The  data  in  Figs.  5  and  6  also  indicate  that  the  distribution 
function  for  stresses  computed  from  stereoimaged  strains  may  not  be  the 
same  as  that  derived  from  subgrain  data.  This  situation  is  still  being 
explored. 
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Figure  lb.  Mohr's  circles  of  stress ,  fine  grid,  dry  nitrogen  environment, 
AK  *  11  MN/rn3/2. 
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Figure  2a.  Mohr's  circles  of  stress,  course  grid,  wet  air  environment, 
AK  =  11  MN/m3/2. 
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Figure  3b.  Mohr's  circles  of  strain,  fine  grid,  dry  nitrogen  environment, 
&K  -  11  MN/m3/2. 
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Figure  4b.  Mohr's  circles  of  strain,  fine  grid,  wet  air  environment, 
AK  ■  11  MN/m3/2. 
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Figure  5.  Effective  stress  vs.  distance  from  crack  tip  for  stress 
derived  from  subgrains  and  from  stereo imaging . 
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Table  1 

A  a 

—  at 
®y 

crack  tip 

Environment 

Subgrain 

Correlation 

Computed  from 
Strains 

Dry  Nitrogen 

2.3 

2.46 

Wet  Air 

2.0 

2.17 

Act 

wet  i 

» __  r*0 

A  a,  1 

dry 

0.89 

0.88 

Calculation  of  mean  stress  range  from  the  principal  stress 
ranges  is  easily  accomplished  using  the  data  of  Figures  1  and  2  for  the 
plane  stress  case.  Figure  7  shows  the  computed  mean  stress  ranges  for  the 
two  environments  directly  ahead  of  the  crack.  The  significance  of  the 
oscillations  is  not  clear,  and  it  is  difficult  to  determine  just  what 
factors  may  be  causing  them  because  of  the  complex  computational  procedure 
used  in  their  calculation.  For  the  "plane  strain"  assumptions,  the  compu¬ 
tational  procedure  is  much  simpler.  Mean  stress  is  computed  directly  from 
the  average  strain  (see  Appendix  B) ;  therefore,  the  average  strain  is  also 
plotted  on  the  figure.  Mean  stress  as  determined  from  average  strain 
depends  on  the  algebraic  sum  of  (ATi  +  AT2>  ,  but  since  AT2  is  nearly  equal 
to  the  negative  of  AT^,  it  is  the  difference  of  two  nearly  equal  numerical 
values  which  fixes  the  value  of  mean  stress  range.  Thus,  small  errors  in 
the  determination  of  AI^  may  be  responsible  for  the  oscillations  observed. 
We  have  yet  to  resolve  whether  or  not  the  oscillations  are  real  or  due  to 
Inaccuracies  in  determining  of  ATi  and  AT2>  At  this  time,  it  appears  as 
though  the  magnitude  of  om  is  relatively  insigificant  as  a  factor  in 
environmentally  assisted  fatigue  crack  growth  in  this  material  except 
perhaps,  very  near  (r  <5  pm)  the  crack  tip. 


In  order  to  obtain  om,  the  maximum  mean  stress,  calculated 
values  of  the  mean  stress  range  Aom  must  be  added  to  omin>  However,  there 
is  good  reason  to  believe  (see  Section  IV)  that  om^n  <<  Aam,  so  that 
A®m  3  ®m  • 


Mean  Stress 


Mean  Strain  (Z) 


III.  CALCULATION  OF  THE  ENERGY  EXPENDED  IN  CREATING  A  UNIT  AREA  OF  NEW 
CRACKED  SURFACE 

A  computational  procedure  has  been  developed  to  make  an  estimate  of 
the  energy  expended  in  creating  new  cracked  surface  from  the  data  of  a 
single  load-unload  pair  of  photographs  of  the  crack  tip.  This  section 
outlines  the  results  of  these  computations  and  compares  them  to  those  pre¬ 
viously  obtained  using  subgrains(i)  and  strain  gages (3)  for  strain  and 
energy  measurement.  A  more  detailed  account  of  the  computational  procedure 
and  accompanying  assumptions  is  given  in  Appendix  C. 

The  energy  computation  begins  with  a  knowledge  of  the  effective 
strain  increment  Aeeff  at  any  point  (x,y)  within  the  crack  tip  plastic 
zone.  This  strain  is  computed  through  Eq.  A8  in  Appendix  A,  from  the 
principal  strains;  these  are  computed  from  ex,  £y  and  eXy  as  determined 
by  stereoimaging  analysis. 

By  plotting  stress  range  vs.  strain  range  for  uncracked  specimens, 
the  cyclic  stress-strain  and  hysteretic  energy  loss  functions  may  be 
obtained.  The  work  per  cycle  Wc  is  computed  as  the  area  within  a  hyster¬ 
esis  loop.  Graphical  representation  of  the  results  are  given  in  Figure  8. 
The  line  through  the  data  fits  the  equation 

W  “  W  'Ae®  (2) 

c  co 

where  Wco'  ■  675  MJ/m^  and  m  “  1.14. 

An  equation  of  the  same  form  as  (2)  may  be  derived  Independently 
using  subgrain  data. (2)  In  this  case  correlation  of  the  work  per  cycle 
with  subgrain  size  gives 

W  -  W  d"1  (3) 

C  CO 

and  correlation  of  Ac  with  subgrain  size  gives 

Ae  -  E  d"m  (4) 

o 

where  m  ■  0.84.  Combining  (3)  and  (4)  yields 

wc  -  wco<r>5  (5) 

o 


where  m  *  1.19.  Comparison  of  the  exponential  values  shows  very  good 
agreement  between  the  directly  measured  and  derived  values.  Using  Wco  ■ 
3x10^  J/m^CD  and  E0  ■  .12,(25  gives 
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The  pre-exponential  factor  as  derived  from  subgrains  is  a  factor  of  about 
20  too  small,  but  the  exponent  is  in  agreement.  It  is  thus  concluded  that 
Eq.  (2)  has  the  correct  functional  form,  and  it  is  this  relation  which  was 
used  together  with  measured  values  of  Aegff  to  compute  the  energy  expended 
in  creating  new  cracked  surface,  *^/dA. 

For  the  stereoimaging  strain  data  previously  presented, (2)  t^e 
computational  results  using  Eq.  (2)  and  the  procedure  of  Appendix  C  are 
given  in  the  following  Table  2. 


Table  2 


Environment 

AK 

MN/m3/2 

Wj 

MJ/m2 

Dry  N2 

11 

11.4 

Wet  Air 

11 

4.2 

7.2  = 

energy 
of  the 

equivalent 

environment 

These  results  are  compared  with  those  previously  obtained  by  use  of  sub¬ 
grains  in  Table  3,  where  an  additional  amount  of  energy  of  3.6xl06  j/m2 
has  been  added  in  order  to  account  for  energy  dissipated  outside  of  the 
subgrain  forming  region,  as  determined  by  Liaw,  et  al.(^) 


Table  3 


Environment 
Dry  N2 
Wet  Air 


AK 

MN/m3/ 2 

Wj 

MJ/m2 

11 

4.7 

11 

3.9 

0.8  *  energy  equivalent 
of  the  environment 


The  magnitudes  compare  well  for  the  wet  air  environment,  and  differ 
somewhat  for  the  dry  nitrogen  environment,  although  the  difference  is  not 
considered  significant,  considering  the  differences  in  the  measuring  tech¬ 
niques.  The  values  as  determined  from  subgrain  and  strain  gage  data  are 
considered  the  more  accurate  because  of  the  time-averaging  of  the  data 
collection  procedures  for  these  calculations;  the  energy  calculation  from 
stereoimaging  is  only  for  one  instance  (cycle)  and  is  not  time  averaged. 
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It  is,  therefore,  only  an  estimate  of  the  average  value.  The  value  of  the 
computation  is  that  it  can  be  obtained  very  easily  in  the  process  of  com¬ 
puting  of  strain  and  stress,  and,  at  least  in  this  case,  gives  a  reasonable 
estimate  of  the  energy  magnitude,  as  compared  with  more  laborious  methods. 


IV.  STRAIN-RATE  EFFECTS 

In  the  1979  Annual  Report,  displacement  diagrams  for  the  unloaded 
crack  versus  no  crack  were  shown  (Figs.  16  and  17).  An  analysis  of  this 
data  revealed  that  compressive  displacements,  relative  to  the  crack  plane, 
existed  in  the  crack  tip  region  for  the  dry  environment,  and  compressive 
displacements  of  about  half  this  magnitude  existed  for  the  wet  environment. 
In  fact,  these  residual  displacements  were  so  large  that  full  loading  of 
the  crack  tip  was  insufficient  to  overcome  them. 

Further  investigation  of  these  observations  has  been  made.  The 
original  observation  concerned  a  crack  tip  loaded  cyclically  at  0.03  Hz. 
Subsequent  investigation  at  0.3,  0.1,  1.0,  and  5  Hz  has  not  indicated 
similar  crack  tip  displacements,  while  duplication  of  the  0.03  Hz  test 
showed  findings  which  were  qualitatively  the  same  as  previously  reported. 

At  the  higher  frequencies  then,  there  is  no  (detectable)  tensile  or 
compressive  residual  displacement  in  the  crack  tip  region,  indicating 
that  under  such  cyclic  conditions,  the  crack  tip  strains  are  close  to 
completely  reversible. 

The  "compressive  displacement"  phenomena  observed  at  0.03  Hz  is 
probably  due  to  cyclic  creep.  Cyclic  creep  effects  caused  by  tensile 
mean  stress  for  plain  carbon  steel  have  been  recently  reported.^)  The 
low  loading  rate  at  0.03  Hz  is  thought  to  give  sufficient  time  for  tensile 
creep  near  the  crack  tip  to  occur.  When  the  crack  is  observed  unloaded, 
the  relative  displacements  of  the  crack  flanks  gives  the  appearance  of 
compression  on  the  crack  plane,  whereas  actually  the  strains  are  tensile. 
These  tensile  strains  are  probably  responsible  for  the  development  of 
compressive  stresses  due  to  the  complete  closure  of  the  crack,  but  these 
have  not  actually  been  computed. 

The  lack  of  development  of  "compressive  strains"  at  the  crack  tip 
for  loading  rates  greater  than  0.03  Hz,  taken  together  with  no  abrupt 
change  in  the  relation  between  crack  growth  rate  and  loading  rate,  indi¬ 
cates  that  this  cyclic  creep  phenomenon  is  not  of  major  significance  in 
the  lengthening  of  the  fatigue  crack,  but,  rather,  this  phenomenon  is 
just  another  manifestation  of  the  strain  rate  dependence  of  plastic 
deformation  of  low-carbon  steel. 
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V.  INTERPRETATION  OF  DATA 

The  information  on  crack  tip  plasticity  which  has  been  developed  on 
this  investigation  to  date  is  reviewed  in  this  section  and  two  interpreta¬ 
tions  of  it  are  made.  The  first  concerns  the  observations  and  measurements 
of  the  energy  dissipated  in  advancing  the  fatigue  crack  and  its  variation 
with  stress  intensity  factor  and  environment.  Part  A  of  this  section  is  a 
reprint  of  a  manuscript  which  was  sent  to,  and  has  been  tentatively 
accepted  by.  Fatigue  of  Engineering  Materials  and  Structures.*  Part  B 
summarizes  the  findings  of  the  investigation  which  are  related  to  the 
detailed  mechanistic  understanding  of  crack  propagation,  and  discusses  the 
possible  mechanisms  of  crack  propagation  which  are  compatible  with  these 
results . 


*Figure  numbers  and  references  for  this  section  are  separate  from  those 
for  the  rest  of  the  report. 
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A.  INCORPORATING  THRESHOLD  AND  ENVIRONMENTAL  EFFECTS  INTO  THE 
DAMAGE  ACCUMULATION  MODEL  FOR  FATIGUE  CRACK  GROWTH 

ABSTRACT 

Measurement  of  the  energy  expended  during  growth  of  a  fatigue  crack 
in  the  near  threshold  region  and  how  it  is  affected  by  a  water  vapor  en¬ 
vironment  is  coupled  with  direct  observation  of  the  crack  tip  to  produce 
a  modification  of  the  damage  accumulation  model  for  fatigue  crack  growth 
which  incorporates  threshold  and  environmental  effects. 

INTRODUCTION 

Although  a  large  number  of  equations  exist  for  calculating  the 
growth  rate  of  a  fatigue  crack,  there  is  no  general  agreement  as  to  which 
formulation  approach  is  the  most  appropriate.  This  situation  is  at  least 
partly  true  because  of  the  very  large  number  of  variables  to  be  con¬ 
sidered  in  any  experimental  or  applications  situation.  Several  good 

recent  catalogues  of  the  different  analytical  approaches  and  the  crack 

(1  7) 

growth  equations  which  result  from  them  have  been  given.  Of  these 

approaches,  the  crack  growth  models  (as  opposed  to  semi-empirical  formu¬ 
lations)  fall  into  essentially  two  categories: 

1)  Those  based  on  a  characteristic  crack  tip  opening 

displacement  (CTOD)  criteria,  which  lead  to  the 

da  2 

general  relationship  —  <*  AK 

aN 

2)  Those  based  on  a  strain  accumulation  or  damage 

criteria,  which  lead  to  the  relationship  —  «  AK 

dN 

where  ^a/dN  is  the  crack  growth  rate,  and  AK  is  the  cyclic  stress  intensity 
factor. 

An  excellent  review  of  the  details  of  derivation  of  the  models 

(3) 

fitting  into  these  categories  is  given  by  Irving  and  McCartney. 

These  authors  point  out  that  values  of  the  exponent  of  AK  have  been 


tffXA  M  - 
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measured  between  2.3  to  19,  but  that  most  reliable  data  taken  in  an  inert 

environment  give  a  value  of  4  +  0.5. 

A  number  of  strain/damage  accumulation  models  have  been  proposed. 
(3) 

Weertman  seems  to  be  the  first  to  have  formulated  such  a  model,  fol¬ 
lowed  shortly  thereafter  by  Rice.  The  theoretical  concept  was  then 
further  developed  by  Weertman^  and  Mura  &  Lin,^  giving  the  general 
relation  for  crack  growth  rate  as: 


da  =  (AK)4 

dN  2 
per  y 


(1) 


where  u  =  shear  modulus,  a =  an  appropriate  measure  of  strength,  F  = 

a  geometric  factor,  and  y  =  the  mechanical  work  input  required  for 

creation  of  a  unit  area  of  new  crack  surface.  Another  way  of  describing 

y  would  be  as  the  effective  surface  energy.  Raju^  has  shown  that  Eq. 

(1)  is  a  special  case  of  the  generalized  relation  he  derived  using  an 

energy  balance  criteria  for  fracture. 

It  should  be  noted  that  Eq.  (1)  was  derived  for  steady-state  crack 

growth  with  the  value  of  AK  constant  or  ever  increasing  and  with  the 

linear,  or  Stage  II,  portion  of  the  crack  growth  (cia/dN  vs  AK)  curve,  and 

(2) 

without  consideration  of  any  environmental  effect.  Irving  and  McCartney 
have  further  refined  Eq.  (1)  to  account  for  the  non-linear  region  of 
^a/dN  vs  AK  as  K  -*  K^  (the  fracture  toughness  stress  intensity) ,  and  to 
account  for  the  effects  of  changing  R  (Kmin/Kmax) .  Very  recently 
Chakrabortty^  has  formulated  a  model  which  relates  Stage  II  fatigue 
crack  growth  to  low  cycle  fatigue  properties,  omitting  the  effects  of 
environment,  and  including  the  threshold  region  only  circumspectly. 
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The  purpose  of  this  paper  is  to  show  how  Eq.  (1)  may  be  interpreted 

in  the  region  of  low  AK,  as  AK  ■+•  AK  (the  threshold  value  of  cyclic  stress 

Trl 

intensity),  and  how  the  effects  of  environment  may  be  included  in  this 
description  of  fatigue  crack  propagation. 


MEASUREMENT  OF  THE  SURFACE  ENERGY 


A  considerable  body  of  experimental  data  for  numerous  materials 
exists  to  support  a  relation  like  that  of  Eq.  (1),  so  it  is  considered  as 
valid  for  steady-state  growth,  but  the  definitions  of  o,  F  and  y  all 
require  some  interpretation.  The  appropriate  definition  of  a  is  probably 
a  "cyclic  yield  stress"  or  something  akin  to  it,  and  F  is  probably  best 
described  (in  practical  terms)  as  a  factor  for  correlating  experimental 
data. 

The  work  required  for  creation  of  a  unit  area  of  crack,  y,  has  been 

measured  for  several  steels  and  aluminum  alloys  by  a  number  of  workers  at 

(8) 

Northwestern  University;  see  the  recent  work  of  Fine''  for  a  compendium 

of  results.  These  measurements  utilize  a  strain  gage  technique,  which 

measures  work  absorption  by  the  specimen  as  the  fatigue  crack  lengthens. 

Due  to  the  size  of  the  strain  gages,  these  measurements  cannot  be  made 

within  200  pm  of  the  crack  plane. 

(9) 

Davidson  &  Lankford  have  made  similar  measurements  on  low-carbon 
steel  using  electron  channeling  contrast  to  observe  the  subgrain  forming 
region  near  the  crack  tip.  A  comparison  of  the  two  techniques  by  Liaw, 

Fine  &  Davidson  determined  that  the  extrapolation  technique  within 
200  pm  of  the  crack  tip,  used  in  applying  the  strain  gage  method,  was  in 
agreement  with  the  measurements  made  in  that  region  by  electron  channeling. 


W 
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and  that  the  values  of  y  derived  by  channeling  were  too  low  (by  about  a 

factor  of  2)  because  considerable  hysteretic  loss  occurs  in  the  specimen 

outside  of  the  subgrain  forming  region. 

The  effect  of  environment  on  the  magnitude  of  y  has  been  determined 

experimentally  by  Davidson  &  Lankford^’ to  be  AK  dependent  in  the  near 

threshold  region,  both  for  inert  and  water  vapor  environments.  Fig.  1. 

The  crack  growth  behavior  observed  fits  the  definition  of  "true  corrosion 

(12)  (13) 

fatigue"  as  defined  by  McEvily  and  Wei,  and  Austen  and  Walker. 

The  data  in  Fig.  1,  together  with  other  experimental  observations,  will 
now  be  examined  and  a  model  presented  for  explaining  the  observed  results. 


A  MODEL  FOR  THE  NEAR-THRESHOLD  REGION 


The  observed  energy  required  for  crack  propagation,  shown  in  Fig.  1, 
gives  a  value  of  y  for  Eq.  (1)  that  is  constant  in  the  limit  (at  higher 
AK) .  The  AK  corresponding  to  this  limiting  value  of  y  constitutes  the 
lower  limit  of  fully  developed  Stage  II  crack  propagation.  Most  of  the 
region  shown  in  Fig.  1  covers  that  part  of  the  crack  growth  curve  where 
the  threshold  AK  for  crack  propagation  is  being  approached. 

The  crack  tip  region  in  this  low-carbon  steel  has  been  observed  in 

(14) 

the  scanning  electron  microscope  (SEM)  using  a  cyclic  loading  stage, 
both  after  fatigue  crack  propagation  in  vacuum  and  in  a  100%  relative 
humidity  environment.  Although  the  SEM  requires  a  vacuum  for  operation, 
there  is  a  residual  effect  of  environment  which  may  last  up  to  a  hundred 
cycles.  Thus,  high  resolution  photographs  of  the  crack  tip  region  may  be 
made  for  both  environments  in  the  SEM.  Data  are  then  derived  from  these 
photographs  by  the  stereoimaging  technique,  which  allows  displacements 


L 

L 
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of  material  in  the  plastic  zone  to  be  both  seen  and  measured,  as  well  as 
values  of  CTOD  and  opening  load  to  be  measured.  The  displacements  may 
also  be  used  to  obtain  five  of  the  nine  elements  of  the  non-symmetric 
strain  tensor. 

Use  of  these  techniques  has  shown  that  the  environment  affects  the 
mode  of  crack  tip  opening,  the  CTOD,  the  opening  load,  and  the  magnitude 
and  distribution  of  strains  in  the  plastic  zone.  Changes  in  stress 

intensity  are  also  observed  to  affect  the  opening  mode. 

The  model  being  proposed  here  links  these  mode  change  observations 
with  the  energy  measurements  given  in  Fig.  1.  The  low  AK  region  is  ex¬ 
panded  in  Fig.  2,  and  the  trends  are  extrapolated  to  encompass  the 
threshold  part  of  the  crack  propagation  curve,  and  to  show  the  elements 
of  the  model.  The  first  observation  is  that  a  limiting  value  of  y  is 
reached  for  AK  in  the  region  where  Eq.  (1)  would  be  expected  to  be  appli¬ 
cable  (and  for  which  it  was  derived). 

For  the  dry  nitrogen  environment,  the  curve  then  develops  a  "U" 
shape  as  AK  is  decreased.  The  higher  AK  portion  of  the  "U"  is  a  transi¬ 
tion  region  where  the  mode  of  crack  opening  is  changing  from  predominantly 
Mode  I.  The  lower  AK  region  is  postulated  to  exhibit  a  mode  of  crack 
opening  which  is  dominated  by  Mode  II.  Since  the  crack  is  driven  forward 
mainly  by  Mode  I  opening,  the  rate  of  crack  growth  in  this  region  is 
decreasing;  i.e.,  the  exponent  for  AK  (as  in  the  Paris  equation  ^  -  CAKm) 
is  becoming  greater  than  4.  As  the  degree  of  Mode  I  decreases,  the  crack 
is  beginning  to  wander  through  the  material,  creating  a  very  rough 
fracture  surface.  Microscopically,  this  motion  of  the  crack  on  planes 
not  perpendicular  to  the  load  axis  results  in  an  increasing  component  of 
Mode  II  opening  as  AK  decreases.  In  the  limit,  with  decreasing  AK,  an 
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ever  increasing  input  of  mechanical  energy  to  the  specimen  is  required 
for  each  unit  area  of  crack  created,  which  is  the  definition  of  a 
threshold  stress  intensity  factor.  This  near  threshold  region  can  be 
viewed  in  another  way:  since  Mode  II  dominates  the  response  of  the  crack 
to  loading,  the  volume  of  material  absorbing  energy  is  small*  and  con¬ 
centrated  approximately  along  the  plane  of  crack  propagation;  the  crack 
is  essentially  made  up  of  short  sections  of  cracks  growing  in  Stage  I, 
a  condition  noted  by  numerous  other  Investigators.  As  AK  is  increased, 
it  is  increasingly  difficult  to  maintain  this  growth  mode,  and  finally  a 

transition  to  Stage  II  begins.  The  details  of  this  transition  depend 

Kmin 

upon  a  number  of  factors:  surely  the  magnitude  of  R  **  'Kmax  grain 

size  will  be  important,  as  perhaps  are  several  other  metallurgical 

variables.  The  postulated  effect  of  increasing  R  is  shown  on  Fig.  2, 

and  it  is  imagined  that  decreasing  the  grain  size  should  have  a  similar 

effect.  Some  experimental  verification  of  the  grain  size  and  R  effects 
(18) 

exists,  although  detailed  analysis  of  these  experimental  findings  in¬ 
dicates  considerable  uncertainty  in  their  interpretation. 

The  above  discussion  has  been  directed  at  the  case  of  crack  growth 
in  an  inert  environment.  For  water  vapor,  an  aggressive  environment,  the 
mechanism  change  on  going  from  the  fully  Stage  II  region  (limit)  down  to 
the  AK  threshold  will  essentially  be  the  same,  except  that  there  is  an 
equivalent  input  of  energy  by  the  environment  which  is  AK  dependent.  The 


Actually  the  whole  specimen  is  absorbing  energy,  hysteretically ,  in 
very  small  amounts  per  cycle.  At  the  crack  tip  a  relatively  larger 
amount  of  energy  is  absorbed  each  cycle,  and  it  is  this  energy  to 


1 

\ 


which  this  statement  refers. 
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transition  region  between  Mode  I  dominated  crack  growth  and  Mode  II 
dominated  growth  will  be  lowered,  as  postulated  on  Fig.  2.  It  must  be 
emphasized  that  in  a  water  vapor  environment  for  the  region  AK  »  8  to  12 
MN/m^2,  experimental  observation has  shown  that  Mode  I  is  the  pre¬ 
dominant  opening  mode,  and  crack  wandering  is  largely  suppressed.  A 

3/2 

departure  from  this  condition  is  beginning  at  AK  =  10  MN/m  ,  so  that 
3/2 

below  AK  -  4-6  MN/m  crack  propagation  might  be  expected  to  be  classi¬ 
fied  as  Stage  I.  The  effect  of  increasing  R  would  again  be,  as  in  the 
case  for  an  inert  environment,  to  extend  to  lower  values  of  AK  the 
transition  to  Stage  I  type  growth.  Since  the  mechanically  input  work  for 
growing  a  crack  in  the  water  vapor  environment  is  less,  for  a  given  AK, 
than  in  the  inert  environment,  crack  propagation  rate  should  be  greater, 
as  is  experimentally  observed,  and  the  ^KThresi10^<j  should  be  lowered  for 
the  water  vapor  environment,  as  experimentally  observed. ^12^  Equation  (1) 
may  be  used  to  incorporate  these  observed  changes  through  the  modification 
of  y. 

As  derived,  Eq.  (1)  incorporates  y  as  a  term  which  is  related  solely 
to  the  amount  of  energy  which  was  absorbed  by  the  body  as  the  crack  was 
being  lengthened.  If  the  equivalent  energy  which  is  supplied  by  the  en¬ 
vironment,  designated  Yg,  is  subtracted  from  y  then  the  effective  energy 

y-Y  causes  the  crack  to  advance  faster  for  a  fixed  AK.  No  method  has 
e 

yet  been  devised  for  independently  measuring  yg. 


Note:  These  extrapolations  do  not  consider  changes  in  the  type  of 

(12) 

corrosion  assisted  fatigue,  as  may  occur  in  some  materials  under 
(18) 


some  conditions. 


The  mechanism  by  which  y  is  altered  by  environment  is  not  known, 

3/2 

but  in  a  high  humidity  environment  at  low  AK(<10  MN/m  ) ,  flat  facets 

appear  on  the  fracture  surfaces.  It  is  likely  that  these  facets 

are  indicative  of  intergranular  fracture.  However,  near  AK^h 

(18} 

fracture  facets  cease  to  appear.  These  observations  taken  together 

with  the  decreased  crack  tip  strain  caused  by  a  water  vapor  environment ^ ^ 
are  consistent  with  a  lowering  of  y. 

Incorporating  the  effects  of  microstructure  on  the  Mode  I/Mode  II 
mix  of  crack  tip  opening  into  Equation  (1)  is  a  more  complex  process,  and 
one  for  which  few  guidelines  exist.  The  lower  the  stacking  fault  energy 
(SFE) ,  the  higher  the  AK  which  might  be  expected  to  be  predominated  by 
Mode  II,  thus  y  could  be  expected  to  increase  with  increasing  SFE. 
Decreasing  the  grain  size  (when  this  is  a  relevant  parameter)  could 
reasonably  be  expected  to  have  about  the  same  effect  as  increasing  R,  and 
the  evidence.  Fig.  3,  shows  this  trend  may  be  correct.  The  effect  of 
other  metallurgical  factors  on  y,  such  as  precipitate  and  dispersoid 
spacing,  texturing,  etc. ,  is  not  known. 

DISCUSSION 

The  theoretically  derived  equation  for  crack  growth,  Eq.  (1),  has 

been  shown  to  have  validity,  at  least  for  some  materials  and  conditions, 

when  steady-state  Stage  II  fatigue  crack  propagation  conditions  are 

(2) 

applicable.  Irving  and  McCartney  have  modified  the  equation  to 

(8) 

accommodate  the  high  AK  region  where  K  -*■  K^,  and  Fine v  has  suggested 
that  y  is  dependent  on  AK,  thereby  accommodating  those  materials  where  the 
exponent  of  AK  <  A.  The  data  of  Figure  1  makes  it  possible  to  modify 
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Eq.  (1)  for  the  condition  where  AK  -*■  T*ie  basis  for  Equation  1 

to  provide  a  framework  for  rationalizing  fatigue  crack  growth  in  numerous 
alloy  systems  seems  to  be  established,  at  least  on  a  conceptual  basis. 
Unfortunately  it  is  difficult  to  measure  y,  and  it  needs  to  be  determined 
for  a  larger  number  of  alloy  systems  under  known  conditions  before  knowing 
the  breadth  of  applicability  of  this  framework. 

Likewise  a  direct  method  for  measurement  of  y  „  is  not 

environment 

available,  and  it  has  been  measured  indirectly  only  for  low-carbon  steel. 

To  make  this  situation  more  complex,  the  crack  opening  mode  is  changing 
in  the  near-threshold  region,  and  this  interacts  with  both  the  environment 
and  the  microstructure  of  the  material.  Because  of  these  complex  inter¬ 
actions  which  will  change  with  alloy,  microstructure  and  environment,  it 
may  be  that  an  approach  which  reflects  these  complexities  through  y  pro¬ 
vides  the  best  overall  description  of  the  threshold  phenomenon.  There  is 
also  hope  that  y  will  become  easier  to  measure  for  a  broader  range  of 
materials  and  environments  as  methods  are  devised  to  use  the  strain  in¬ 
formation  derived  by  the  stereoimaging  method  for  calculating  the  work 
required  for  the  creation  of  new  fracture  surface. 

CONCLUSIONS 


1)  The  effect  of  decreasing  stress  intensity  factor  magnitude  on 
the  mechanical  work  expended  in  forming  a  unit  area  of  new 
fatigue  crack  surface  may  be  explained  by  changes  in  the  pre¬ 
dominant  mode  of  crack  opening  from  Mode  I  to  Mode  II. 


2)  The  value  of  ^pyiresh0id  by  postulating  that  at 

very  low  values  of  AK,  the  change  to  Mode  II  is  nearly  com¬ 
plete  and  that  the  magnitude  of  loading  is  so  small  that  the 
large  amount  of  work  which  must  be  put  into  the  specimen 
results  in  a  vanishingly  small  crack  growth. 

3)  The  effect  of  environment,  for  a  decreasing  AK,  is  to  move 
the  transition  of  crack  opening  mode  from  Mode  I  to  Mode  II 
down  to  lower  values  of  AK. 

4)  When  the  above  observations  are  brought  together  and  expressed 

in  the  form  of  the  energy  required  to  make  a  unit  area  of  new 

cracked  surface,  they  fit  conveniently  into  the  crack  growth 

(3) 

rate  equation  originally  formulated  by  Weertman  for  Stage  II 
fatigue  crack  propagation.  With  the  modification  of  Irving 
and  McCartney  for  the  high  AK  region  this  type  of  equation  may 
be  used  over  the  full  range  of  crack  growth  rates  as  a  frame¬ 
work  for  calculating  steady-state  fatigue  crack  growth. 

5)  The  expected  effect  of  changing  R  ratio  is  shown  in  Fig.  2, 

but  the  change  in  AKT^re8h0ld  with  changes  in  both  R  ratio 

(18) 

and  grain  size  is  likely  to  be  very  complex. 
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Figure  1 


Energy  expended  Co  grow  a  fatigue  crack  as  a  function  of 
cyclic  stress  intensity  showing  the  effect  of  environment 


(EFFECTIVE  SURFACE  ENERGY)  J/m 


The  effect  of  grain  size  and  R  ratio  on  threshold  AK  for 
16  different  steels  tested  in  air.  The  lines  describe  the 
bounds  for  the  correlation,  which  includes  102  data  points. 
After  S.  G.  Druce. <18> 
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B.  Summary  and  Interpretation  of  Information  Related  to  the 

Mechanism  of  Environmentally  Assisted  Fatigue  Crack  Propagation 

To  date,  experimental  observations  have  determined  that  the 
effects  of  a  water  vapor  environment  on  crack  tip  plasticity  parameters 
for  low-carbon  steel  are  as  follows: 


1) 

The  energy  dissipated  in  forming  new  crack  surface  is 
lowered,  and  is  AK  dependent  in  the  near  threshold 
region. 

2) 

The  stress  range  magnitude  is  lowered,  and  is  AK 
dependent. 

3) 

The  strain  range  magnitude  is  lowered. 

4) 

The  size  of  the  subgrain  forming  region,  which  is 
proportional  to  the  plastic  zone  size,  is  decreased, 
and  is  AK  dependent. 

5) 

The  mix  of  deformation  modes  is  changed  to  more  Mode  I 
(less  Mode  II) . 

6) 

On  the  fracture  surface,  "brittle,"  probably  inter¬ 
granular  facets,  begin  to  appear  as  AK  is  decreased. 

From  these  observed  changes,  it  is  clear  that  both  the  magnitude  and 
spatial  extent  of  crack  tip  plasticity  are  reduced  by  a  water  vapor 
environment  as  AKth  is  approached.  But,  what  do  these  findings  imply 
about  the  material  property  changes  in  the  element  which  will  be  fractured 
thereby  advancing  the  crack,  on  the  next  tensile  portion  of  the  loading 
cycle?  And,  exactly  where  are  the  material  properties  affected  by  the 
environment — It  seems  clear  enough  that  they  are  affected  just  adjacent 
to  the  crack  tip,  but  are  they  also  affected  farther  out  in  the  plastic 

zone? 


The  only  atomic  species  available  at  the  crack  tip  in  a  water 
vapor  environment  are  hydrogen  and  oxygen.  From  the  work  of  Dwyer, 
Simmons  and  Wei, (6)  there  seems  to  be  little  doubt  that  hydrogen  is  the 
element  responsible  for  whatever  changes  are  found.  Hydrogen  could  alter 
material  properties  in  one  or  more  of  four  ways.  It  could: 

1)  increase  or  decrease  the  flow  stress; 

2)  increase  or  decrease  the  strain  to  fracture; 

3)  change  the  work-hardening  characteristics;  and 

4)  lower  the  surface  energy  at  the  crack  tip. 

For  material  properties  to  be  altered  beyond  the  immediate  crack  tip 
region,  farther  out  into  the  plastic  zone,  transport  of  the  hydrogen  into 
that  region  would  be  required.  This  movement  of  the  hydrogen  could. 


conceptually,  be  accomplished  by  concentration  gradient  and  stress  gradient 
driven  diffusion  and  by  dislocation  sweep  in. 

First,  the  possible  material  property  changes  immediately 
adjacent  to  the  crack  tip  will  be  examined: 

1)  Since  the  energy  dissipated  in  forming  the  crack 
is  lowered,  the  area  beneath  the  stress-strain 
curve  must  be  decreased. 

2)  Since  the  crack  tip  strain  magnitude  is  decreased 
by  about  a  factor  of  2,  this  must  be  a  major  con¬ 
tributor  to  the  decrease  in  energy  dissipation 
(and  area  beneath  the  stress-strain  curve) . 

3)  The  decrease  in  stress  magnitude  determined  both 
from  subgrain  size  and  from  strains  by  stereo¬ 
imaging  assumes  no  change  in  the  stress-subgrain 
or  stress-strain  relations  for  hydrogen  charged 
material.  If  this  assumption  was  not  correct, 
then,  effectively,  the  correlation  between  stress 
range  and  subgrain  size  would  have  to  be  thought 
of  as  being  altered  by  the  presence  of  hydrogen. 

Indeed,  it  is  possible  to  make  the  assumption  that 
the  stress  range  at  the  crack  tip  is  not  altered 
by  environment,  and  back  calculate  the  functional 
relation  between  stress  range  and  subgrain  size. 

But  since  the  alteration  of  stress  range  at  the 
crack  tip  is  AK  dependent, this  would  mean 
that  a  different  correlation  between  stress  range 
and  subgrain  size  would  be  required  for  each  AK, 
and  this  is  not  physically  plausible.  Thus,  the 
decrease  in  the  stress  range  in  a  wet  environment 
appears  to  be  a  real  effect,  and  must  be  explained 
by  either  a  change  in  yield  stress  or  work 
hardening  or  both.  If  the  water  vapor  environment 
is  not  to  alter  the  ratio  of  stress  range/yield 
stress,  then  the  effect  of  hydrogen  must  be  to 
decrease  the  yield  stress,  which  is  not  likely  for 
low-carbon  steel^>®)  (although  it  might  be  for 
pure  iron) .  A  decrease  in  yield  stress  should 
cause  a  larger  plastic  zone,  but  a  smaller  plastic 
zone  is  actually  observed,  indicating  that  the 
yield  stress  ought  to  be  increasing.  If  yield 
stress  increases,  then  this  would  cause  the  stress 
range/yield  stress  ratio  to  be  further  lowered. 

4)  The  possibility  of  a  change  in  work  hardening 
coefficient  may  also  be  considered  by  examining 
the  stress  and  strain  distributions.  A  change  in 
work  hardening  would  be  manifested  in  a  change  in 
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stress  and  strain  distribution  functions  in  the 
way  outlined  by  Rice  and  Rosengren, (9)  who  found 
^  1  / 1  \ 
a  «  r  i+N  and  e  «  r  i+n.  Previously,'-1''  the 
stress  distribution  was  found  to  be  a  «  r~-26; 
this  value  came  from  the  correlation  of  stress 
range  and  subgrain  size.  From  this,  a  strain 
distribution  e«  r-< was  derived;  it  was  possible 
to  obtain  a  check  on  the  validity  of  this  by 
stereoimaging. (2)  These  distribution  functions 
would  result  if  N  ■  0.31,  which  is  about  the  value 
obtained  by  Landgraf^  for  low-carbon  steel.  The 
effect  of  water  vapor  is  to  change  the  magnitudes 
of  these  functions  at  the  crack  tip,  but  not  the 
distribution  function  (slope) .  It  is  therefore 
concluded  that  work  hardening  of  the  material  is 
not  changed  substantially  by  the  presence  of 
hydrogen. 


5)  To  examine  the  possibility  that  hydrogen  lowers 

the  surface  energy  at  the  crack  tip  and  that  this 
is  the  relevant  material  alteration,  it  is  useful 
to  compare  the  magnitudes  of  the  energy  dissipated 
in  lengthening  the  crack  versus  the  surface  energy. 
Table  4.(2) 


Table  4 


_ Condition _ 

No  Environmental  Effect 
With  Environment 
Surface  Energy 


Energy  (J/m2) 
17.1xl05 
1.5x10s 
2 


AK(MN/m3/2) 

8 

8 


Thus,  the  environmental  effect  is  in  the  range  of 
15. 6xl03  J/m2  as  compared  to  2  J/m2  for  the  surface 
energy.  If,  for  example,  the  environment  were  to 
lower  the  surface  energy  to  1  J/m2,  no  model  exists 
which  can  explain  how  this  might  have  a  1x10^  J/m2 
effect  on  the  energy  dissipated  in  forming  new 
crack.  This  possibility,  therefore,  is  excluded. 


The  above  analysis  argues  for  the  fact  that  hydrogen  alters 
only  one  material  property  in  a  substantial  way;  the  strain  to  fracture 
is  decreased  by  about  a  factor  of  2.  The  yield  and  flow  stress  may,  or 
may  not,  be  altered,  but  if  there  is  any  change,  it  is  not  very  large; 
not  over  +  15%. 
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The  conclusion  that  work  hardening  is  not  changed  could  also 
be  interpreted  as  follows:  since  the  distribution  functions  of  stress 
and  strain  are  not  substantially  altered,  but  the  magnitude  (coefficient) 
of  the  function  is,  this  is  indicative  of  only  a  small  volume  of  the 
material  immediately  at  the  crack  tip  which  has  been  affected  by  the  hy¬ 
drogen,  rather  than  a  substantial  portion  of  plastic  zone.  If  correct, 
then, this  finding  could  mean  that  the  effect  of  hydrogen  was  simply 
"embrittlement,"  and  would  obviate  the  requirement  for  hydrogen  to  be 
transported  into  the  material  to  any  great  distance  by  either  diffusion 
or  dislocation  sweep-in. 


VI.  SUMMARY  AND  CONCLUSIONS 

1.  A  computational  procedure  for  calculating  stress  from  strain, 
(as  determined  by  stereo imaging)  has  been  found,  and  stresses 
in  the  crack  tip  plastic  zone  have  been  computed. 

2.  The  magnitude  of  the  crack  tip  stresses  calculated  from 
measured  strains  agrees  reasonably  well  with  those  determined 
previously  using  subgrains,  as  detected  by  electron  channeling. 
The  distribution  of  stress  along  some  radial  directions,  how¬ 
ever,  differs  between  the  two  methods. 

3.  The  strains  determined  by  stereo imaging  have  also  been  used 
to  calculate  the  energy  dissipated  in  forming  new  crack 
surface,  and  the  results  agree  reasonably  well  with  those 
obtained  from  subgrain  size  distribution,  as  previously 
determined . 

4.  Previously  reported  crack  tip  "clamping  displacements"  have 
been  further  investigated  and  found  to  be  loading  rate 
dependent.  As  such,  they  are  thought  to  be  a  cyclic  creep 
effect,  and  not  directly  related  to  fatigue  crack  propagation. 

5.  The  previously  determined  changes  of  energy  dissipation  to 
form  new  crack  surface  with  environment  and  AK  have  been 
incorporated  into  the  damage  accumulation  model  for  fatigue 
crack  growth  and  a  manuscript  prepared  for  publication. 

6.  The  crack  tip  plasticity  data  which  has  been  collected  to 
date  is  summarized,  and  analysed  to  determine  the  effect  of 
the  environment  on  material  properties  in  the  crack  tip 
region.  Hydrogen  is  considered  to  be  the  dominant  environ¬ 
mental  specie.  It  is  concluded,  based  on  the  available 
evidence,  that  the  major  effect  of  hydrogen  is  to  reduce 
the  strain  to  fracture  in  the  region  just  adjacent  to  the 
crack  tip,  and  that  the  material  properties  farther  out  in 
the  plastic  zone  are  not  greatly  affected. 
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APPENDIX  A 

CALCULATION  OF  STRESS 


Calculation  of  the  principal  stress  ranges  from  the  principal  strain 
ranges  has  been  done  according  to  the  method  given  by  Lubahn  and  Felgar.d' 
The  authors  are  endebted  to  Prof.  S.  R.  Bodner  of  the  Technion  for  pointing 
out  this  and  other  references  to  them. 

In  order  to  make  the  calculation,  it  is  necessary  to  make  a  number 
of  assumptions: 

1)  Plane  stress,  03  =  0. 

2)  ri  “  ei  +  ei»  ei  M  plastic  strain,  e^  **  elastic  strain  (Al) 

3)  The  deformation  theory  of  plasticity  pertains,  so  that,  with 
03  -  0, 


el  *  e2 


e2  -  e3 


c3  ~  £1 

-o. 


(A2) 


4)  During  deformation,  plastic  strains  do  not  change  the  volume, 
so  that 


el  +  e2  +  e3 


(A3) 


5)  The  elasticity  relations  are 


<31  y 

ei  =  E  "  E  °2 


<3  2  h 
E  "  E  ° 


■j'VVj 


<A4) 


where  y  *  Poisson's  ratio  and  E  *  Young's  Modulus. 
These  equations  may  be  combined  and  reduced  to  the  quadratic 

2 

a o2  +  ba2  +  c  -  0 

where  a  -  l-2y;  b  ■  -E(2Ti  +  T2);  c  -  -[a^(l-2y)  -  Eoi(2T2  +  T^)] 
so  that 

-b-/b^  -  4ac 
°2  *  2a 


(A5) 


(A6) 
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Measured  values  are  AT^,  Che  principal  sCrain  ranges.  It  is  assumed  chac 
Che  above  equaclons  can  also  be  used  Co  calculaCe  Aa-j_. 

The  compuCaClonal  procedure  for  calculacion  of  principal  sCress 
range  is  as  follows: 

1.  Assume  a  value  of  Aaj_. 

2.  CompuCe  A02  from  (A6) . 

3.  Compute 

s  2  2 

Aoeff  -  (Ac^  -  Aa1‘Aa2  +  Aa2  )  (A7) 

4.  CompuCe  Ae^  from  (A4) . 

5.  CompuCe  Ac^  from  (Al) . 

6 .  CompuCe 

1/2 

Aeeff  "  y^^el2  +  AelAe2  +  (A8) 

7.  DeCermine  Aoeffe  from  che  cyclic  sCress-sCrain  curve* 


m 


(A9) 


8.  Compare  AaeffS  to  Aaeffe.  If  unequal,  assume  another  value  of 
Ao^  and  iteraCe.  A  new  value  of  Ao^  is  determined  from 


A  a. 


eff 


(A10) 


Cyclic  stress-strain  curves  are  usually  obtained  from  a  specimen  stressed 
along  one  axis,  while  che  fatigue  crack  creates  a  biaxial  stress  state. 

It  has  been  assumed  that  the  cyclic  stress-strain  curve  may  be  related  to 
the  multiaxial  case  through  use  of  the  effective  stress  and  strain.  This 
assumption  has  been  examined  by  Lambda  and  Sidebottom, (2,3)  who  found  no 
reason  to  invalidate,  or  greatly  modify,  it. 
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APPENDIX  B 

COMPUTATION  OF  MEAN  (HYDROSTATIC)  STRESS  RANGE 


I. 


2. 


3. 


1 

By  definition  Aom  *  3(Ao^  +  Ao 2> ;  A03  *  0  (Bl) 

where  values  of  Lai  and  La 2  are  determined  from  the  plane 
stress  assumptions  of  Appendix  A. 


If  a  condition  of  plane  (plastic)  strain  is  assumed,  so 
that  AS3  ■  0,  a  value  of  the  mean  stress  can  be  determined 
(but  not  01  and  02  individually).  Ao^  +  Ao£  *  2403. 


(l-2p) 


AT.  +  AT, 
<— *-= - “) 


(B2) 


If  a  condition  of  plane  (total)  strain  is  assumed,  so  that 
AT3  *  0,  Ae3  ■  0,  Ae3  *  0,  then 

Adi  +  Ao2  -  ha3 


1.523Ao3 


where  Ao3  is  given  by  Eq.  (B2) . 

Determination  of  A03  and  La 2  may  be  intractable,  an  extended 
effort  has  not  been  made  to  obtain  solutions  for  them. 


4. 


If  a  condition  of  plane  (total)  strain  is  assumed,  so  that 
AT3  **  0  and  AE3  ■  -Ae3,  then 


E 

<l-2u) 


+  AT, 


) 


Table  Bl 

Summary  of  Plane  Strain  Assumptions 


Designation 

A 

B 

C 


Assumptions 


Ae3  “  0 

L?3  *  0,  Ae^  ■  0,  Ae^ 
ATj  •  0,  Ae3  ■  -Ae^ 


Mean  Stress 


A  a, 


A^  +  AT 


m 


.  £1  /  A 

Aa3  '  (ir^7)( - 2  } 


-  0  1.523AOi 

— Ao 

n 


A  R 

Thus,  the  assumptions  giving  Aom  lie  between  those  yielding  Aom  and 
Ac%c: 


2.  A 
tA  a 
3  m 


.  B  .  3 
A  0  *  -rAo 

m  2  m 


A 


All  these  plane  strain  assumptions  depend  on  the  quantity  (AI*i  +  AT2> , 
but  algebraically,  since  AT2  is  negative,  this  becomes  the  numerical 
difference  between  the  magnitudes  of  AT^  and  AT2,  which  is  usually  a  small 
quantity.  The  term  . . .  »  4.5x10^  MPa,*  thus  the  actual  value  of  Aom  is 


l-2v 


The  mean  strain 


extremely  sensitive  to  small  differences  in  ATi  and  AT2* 

AfL  +  _AT2) 

range  2  is  shown  as  the  number  at  each  point  in  Figures  3  and  4. 

Magnitudes  for  the  plane  strain  estimates  of  Aom  given  in  Table  Bl  may  be 
made  using  these  mean  strain  values  directly.  A  (typical)  mean  strain  of 
0.4Z,  therefore,  gives  Aom  ■  1800+^®®  MPa,  all  of  which  are  believed  to  be 
unrealistically  high. 


Reference 
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^Poisson's  ratio  was  taken  as  y  ■  0.28;  see  Nadi^  for  more  information. 
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APPENDIX  C 

THE  ENERGY  EXPENDED  IN  FORMING  A  UNIT  AREA  OF  NEW  CRACK  SURFACE 


This  appendix  outlines  the  method  used  in  computing  the  energy 
expended  in  creating  a  unit  area  of  new  crack  surface,  Wf,  using  the 
strains  computed  by  stereoimaging  analysis. 


1.  In  the  reduction  of  data  from  a  pair  of  photographs,  one  of 
the  crack  tip  loaded,  the  other  unloaded,  the  strain  ranges 
Aex,  Aey  and  Ayxy  are  determined  at  numerous  points  x,  y 
within  the  plastic  zone  of  the  crack  tip. 


2. 


3. 


From  these  normal  strain  ranges  the  principal  strain  ranges 
are  computed,  using  a  Mohr's  circle  construct.  An  "effective 
strain  range" (1)  is  then  computed  as 


i/: 


aEeff  "  73(Ae1  +  Ae1Ae2  +  Ae2  ) 


(Cl) 


The  hysteretic  energy  loss  per  cycle  is  found  to  be 


W  •  W  'Ae 
c  CO 


m 


Where  WCQ'  -  675  MJ/m^  and  m  »  1.14. 


(C2) 


4.  The  computed  value  of  Aeeff  is  then  put  into  (C2)  and  Wc  is 
obtained  for  each  x,y  point.  A  point  actually  represents  a 
small  volume  of  material,  defined  by  the  x  and  y  spacing  of 
points  on  the  surface  as  A(x,y),  and  one  unit  deep. 

5.  The  total  energy  expended  during  that  one  loading  of  the 
crack  tip  is  the  sum  of  all  the  products,  W£A 


W  -  2  W  (x,y)A(x,y)  (J/m)  (C3) 

x,y  c 

6.  To  obtain  ^/dA,  several  assumptions  are  required: 

a.  For  a  given  value  of  x,  the  function  Aeeff(y)  is 
typical  of  the  strain  history  experienced  by  a 
small  volume  of  material  during  the  passage  of 
the  fatigue  crack. 

b.  The  small  volume  of  material  Ax,  Ay,  Az  represented 
by  each  point  x,y  experiences  a  number  of  cycles 

AN,  which  is  determined  by  Ay/— . 

dN 


T  T*  ’ 
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c.  The  energy  expended  by  each  element  x,y  as  the  crack 
passes  Is  IAN  Aeeff(y)* 

7.  With  these  assumptions:  ^W/dA  ■  I  Z  AN  Ae  ff(x,y).  But, 

x  y  ett 

Z  AN  -  N,  where  N  ■  the  number  of  cycles  which  caused  the 
y 

crack  to  advance  a  distance  from  the  crack  tip  to  the  limit 
of  the  analysis  (which  is  usually  the  edge  of  the  photograph) . 
Therefore,  “w/dA  becomes  simply 

dW/dA  -  (J/m2)  (C4) 

da/dN 

where  W  is  determined  by  Eq.  (C3) . 
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